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ABSTRACT 

In this paper, we develop a full statistical method for the pairwise velocity estimator previously proposed, and apply Cosmicflows- 
2 catalogue to this method to constrain cosmology. We first calculate the covariance matrix for line-of-sight velocities for a given 
catalogue, and then simulate the mock full-sky surveys from it, and then calculate the variance for the pairwise velocity field. By 
applying the 8315 independent galaxy samples and compressed 5224 group samples from Cosmicflows-2 catalogue to this statistical 
method, we find that the joint constraint on and erg is completely consistent with the WMAP 9-year and Planck 2015 best-fitting 
cosmology. Currently, there is no evidence for the modified gravity models or any dynamic dark energy models from this practice, 
and the error-bars need to be reduced in order to provide any concrete evidence against/to support ACDM cosmology. 

Key words. Methods: data analysis, statistical; Cosmology: large-scale structure of Universe, cosmic distance, observations; Galax¬ 
ies: kinematic and dynamics 


1. Introduction 


The study of peculiar velocity field is a powerful tool for prob¬ 
ing the large-scale structures of the Universe. This is because in 
the standard A cold dark matter (ACDM) cosmology, the gravi¬ 
tational instability causes the density perturbations to grow and 
the peculiar velocity field to emerge. On large scales, the pe¬ 
culiar velocity field is directly related to the underlying matter 
perturbations that can be used to test the growth of structure in 
the standard ACDM universe. 

There have been various approaches to using the pecu¬ 
liar velocity field to study cosmology and probe the growth 
of structure. Since in the linear perturbation theory, the veloc¬ 
ity field in real space is related to th e integral of m atter den¬ 
sity contrast with a Newtonian kernel (lPeeblesl|1993h . there has 
been wide interest in comparing the measured velocity field 
with the reconstructed density field from galaxy redshift sur¬ 
veys and testing the linear relation between the two. One of 
the approaches is to reconstruct the linear velocity field from 
the d ensity field and compare it with the measured velocity 
field ( Branchini et alJ l2001t iPike & HudsonI I2005I: iDayis et all 
|201 It Ma et al.ll2012l) . The other approach, the “POTENT”, is 
to use the reyerse process, i.e. reconstructing the grayitational 
potential and density field from the yelo city field and comparing 


potential ana aensity neia trom the yelo city neia ana comparing 
with the measured galaxy density field (iDekel et al.|[l9^ ,ll999l: 
[Hudson et al.l Il99-5t ISigad et all Il998t Branchini et al.l 20001) . 


The results from these practices show that the linear pertur¬ 
bation theory works yery well on scales of 10-100/I 'Mpc, 
and the fitted growth rate factors (/erg) are consistent with the 
ACDM cosmology at low redshifts. The second method is to 
reconstruct the cosmic bulk flow on yarious depths of the lo- 
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cal Uniyerse, which are only sensitiye to cosmological pertur- 


bations on large scales (Kashlinsky et al. 

2008 

Watkins et al.l 

120091 iFeldman et al.ll2010: Watkins & Feldman 

I20L51J). In re- 


cent years, there haye been a few studies that claim to find 
yery large bulk flows on a scale of lOO/i *Mpc or on deeper 
scales that seem to exceed the ACDM prediction by a 3cr con- 


fidence leyel (CL) ( 

Kashlinsky et alJ 20081 IWatkins et alJl2009l: 

IFeldman et alJ201C 

;lMacaulay et al.l201 li 20121). But later stud- 


ies show that this might be due to the systematics that arose 
when co mbining different catalogues with different calibration 
schemes (iNusser & DayisI 1201 ll : iMa & Scott! l2013h . The third 
method is to directly fit the yeloci ty field power spectrum from 
the peculiar yeloci ty field data dMacaulay et alJ 1201 iL 120121 : 
[Johnson et al.ll201^ . The recent results from the six-degree-field 
galaxy suryey data (6dE) show that the fitted yalues of structure 
growth r ate at low redshift are consistent with Planck 2013 cos¬ 
mology dPlanck Collaboration XVlll2014i) . 

In this paper, we consider a different estimator of peculiar ye¬ 
locity field, namely the mean relatiye pairwise yelocity of galax¬ 
ies vi 2 , which is defined as the mean yalue of th e peculiar ye¬ 
locity difference of a galaxy pair at separation r (iFerreira et al.l 
1 19991) . In the fluid limit, th e pairwise yelocity becom es a density- 
weighted relatiye yelocity (lJuszkiewicz et al.ll 19981) . 


Vl2(0 = <Vl - V2>p = 


((vi - V 2 )(l + di)(l + 62)) 

l+/(r) 


( 1 ) 


where v and 6 are the peculiar yelocity and density contrast, re- 
spectiyely, and / is the two-point correlation function. Since the 
line-of- sight yelocities of disc rete galaxies are measured for each 
sample, IFerreira et al.l (1 19991) proposed that the estimator of the 
pairwise yelocity is 


vnir) = [2 Y^(sa - sb)pab] / Pab] > ( 2 ) 
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where - \a ■ i'A is the line-of-sight velocity of galaxy A, 
Pab - r ■ (fA + Fb) is the geometric factor, and the summa¬ 
tion is for all pai rs within a di s tance separation bin Ar. It was 
first proposed in iFerreira et al.l (Il999h that the estimator mea¬ 
sures the cosmological den s ity pa rameter Q, and later it was 
found in lJuszkiewicz et al.l ( |2000|) that the matter density of 
the Universe is close to 0.35 and that the Einstein-de Sitter 
model Q = 1 is inconsistent with the data. The result of a 
low-de nsity universe was further confirmed by iFeldman et al.l 
(l2003h with the Mark III catalogue, Spiral Field I-Band (SFI) 
catalogue. Nearby Early-type Galaxies Survey (ENEAR) cata¬ 
logue, and the Revised Flat Galaxy Catalog (RFGC). This be¬ 
comes the early measurement of matter content of the Uni¬ 
verse before the experiment of the cosmic microwave back¬ 
ground radiation from Wilkinson Microwave Anisotropy Probe 
(WMAP) and strong ly indicates that t h ere is A in the cosmic 
budget. In addition, lJuszkiewicz et alJ (1 1 9991) investigated the 
dynamics of the pairwise motion by calibrating the mean pair¬ 
wise velocity with N-body simulations. They provide a theoret¬ 
ical formula that is very consistent with the N-body si mulation 
result. Additionally, iBhattacharya & Kosowskyl (l2007l) forecast 
the prospectiye constraints on cosmological parameters from the 
pairwise kinetic Sunyaey-Zeldoyich effect (proportional to ye- 
locity fi eld). In 2012, by app lying the pairwise momentum esti¬ 
mator of iFerreira et al.l (119991) into the temperature map, the Ata¬ 
cama Cosmology Telescope team (ACT) proy ided the first detec - 
tion of the kin etic Sunyaey-Zeldoyich effect dHand et alj|2012l) . 
More recently, [Planck Collaboration XXX VUl (l2015h ~haye esti¬ 
mated the pairwise momentum of the kSZ temperature fluctua¬ 
tions of Planck maps at the positions of the Central Galaxy Cata¬ 
logue samples extracted from Sloan Digital Sky Suryey (SDSS- 
DR7) data. They find a ~ 2a CL detection of the kSZ sig- 
n al, which is consis tent and slightly lower than the one found 
in iHand et al.l (l2Q12l) . 

The literature listed aboye about the pairwise yelocity esti¬ 
mator is for the old data set, whereas this work uses a new data 
set, the new compiled Cosmicflows-2 catalogue, to constrain 
cosmology. In addition, we deyeloped a new method of com¬ 
puting the coyariance matrix of the pairwise yelocity field and 
formed a likelihood that directly relates the models of pairwise 
yelocity field with the data. The paper is organized as follows. In 
Sect.|2] we introduce the new Cosmicflows-2 data set. In Sect. [3 
we introduce our statistical method of computing the likelihood 
of pairwise yelocities. In Sect.|4] we present our results of con¬ 
straints and compare them with WMAP nine-year cosmology and 
PlancklQlS cosmology. The conclusion is presented in the last 
section. 

Throughout the paper, unless otherwise stated, 
we use Planck2Ql5 best-fitting cosmological parame¬ 
ters dPlanck Collaboration Xlllj l2015b , i.e. {Qm, erg, h, Qb, 
ln(10'%)} = {0.309,0.816,0.677,0.049, 3.064}. 

2. Data set 
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Fig. 1. Histogram of the line-of-sight peculiar velocity difference in the 
GROUP catalogue between th e calculation by u sing Eq. l[3) (labelled as 
VMa) and the values given in iTullv et all (1201 3h (labelled as Vruiiy)- 

tuation, the tip of the red giant branch (TRGB), or Cepheids. 
The Cosmicflows-2 samples are calibrated at their zero points 
by using two different approaches. One is to use the Cepheid 
period-luminosity relation (Cepheid PER), and the other to use 
the luminosities of red giant branch stars at the onset of core he¬ 
lium burning, at a location in a stellar colour-magnitude diagram 
(CMD) known as the TRGB. 

Eventually, the correlation between the galaxy luminos¬ 
ity and HI line width w as refitted so zero points were deter¬ 
mined (iTullv et al.ll2013h . In addition, iTuIIv & CourtoisI (l2QI2l) 
have demonstrated that no extra systematics were found after 
the zero-point calibration. The results of the two methods of re¬ 
calibrations were compa red, which confir med that the distance 
estimation is unbiased (iTullv et al.l l2013h . We obtained these 
samples from the “VizieR” astronomy data base, which are avail¬ 
able in two tables. Table I provides the entry for every galaxy 
with a distance, consisting of 8315 galaxies in total. Table II 
condenses the galaxies in each group and provides a distance 
for each group, therefore consisting 5224 group entries. In the 
following analysis, we use both tables in the pairwise velocity 
analysis, so we name Table I as GALAXY and Table II as GROUP in 
the following text. 

We made one adjustment of the observed peculiar velocitie s 
in the Cosmicflows-2 samples. On page 20 of iTullv et al.l (1201 3l) . 
it is stated that the line-of-sight peculiar velocity is calculated 
through c z - Hod, where d is the meas ured distance. However, as 
shown in iDavis & Scrimgeouil (1201 4l) and IWatkins & Feldm^ 
(1201 5al) , this estimator makes the velocity estimate biased at the 
level Avp ~ lOOkms^' at z ~ 0.04 and even more biased at 
higher redshifts. The m ore accurate formula for ca lculating the 
line-of-sight velocity is (iDavis & Scrimgeourll2014l) 



The Cosmicflows-2 catalogue (ITullv et al.ll2013l) is a compiled 
catalogue of distances and peculiar velocities of more than 8 000 
galaxy samples. Some of these samples are from new measure¬ 
ments, while others are taken from the literature. The majority 
of the distance s of the samples are m easured through the Tully- 
Fisher re lation dTullv & Fishedll97^ or fundamental-plane (FP) 
relation dPiorgovski & Da^ 19871 i^mpbell et al.ll2014l) with 
roughly 20 per cent of the error for the distances. But there are 
a small portion of the samples whose distances are measured 
from Type-la supernovae light curves, surface brightness fluc- 


where z is the measured redshift, and z is the redshift for the 
unperturbed background. One can obtain z by inverting the mea¬ 
sured luminosity distance 

where H(z) - Hq + z)^ + O-h is the Hubble parame¬ 

ter. In Fig. [1] we show the measured velocity difference of 
GROUP catalogue between Eq. (|3l and the Col. 20 of Table II 
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Fig. 2. Leyi-Histogram of the measured velocity error of GALAXT (Orange) and GROUP (blue) catalogues. The two histograms peak at 500-1 000 
and 1 OOOkms^', respectively. The bin width of the histogram is A(T(Vp) = 200kms“'. f^igAf-Histogram of redshift distribution of GALAXY and 
GROUP samples. The bin width of the plot is cAz = 2 000 km s“'. 



Fig. 3. Full-sky GALAXY (8315 samples) and GROUP (5224 samples) catalogues plotted in Galactic coordinates. The red (light grey) points are 
moving away from us, and the blue (dark grey) points towards us. The size of the points is proportional to the magnitude of the line-of-sight 
peculiar velocity. 


in iTullv et ^ (l2013h . One can see that by adjusting the pecu¬ 
liar velocity calculation, there is a trend towards lower veloc¬ 
ity valueSj_^ich_ii^icates that the original velocities provided 
bv lTullv et alJ (1201 3h are biased towards higher velocity values. 
This is consistent wi t h the pr ediction shown in Figs. 2 and 3 in 
iDavis & Scrimgeou^ (l2Q14tl*[ In addition, we propagate the dis¬ 
tance error into the peculiar velocities and plot the histogram of 
the rms velocity errors in the upper panel of Fig. |2l One can 
see that the two histograms peak at slightly different values, i.e. 
500-1 000 and 1 000 km s“\ respectively. We use the propagated 
velocity errors as the measurement error in the likelihood analy¬ 
sis. 

We plot the redshift distribution of GALAXY and 
GROUP catalogues in the lower panel of Fig. |3 One can see that 
the bulk of the samples are in the range of 0-15 OOOkms 
while a few samples reside at higher redshifts. The two cata¬ 
logues are probing the same redshift range. The GROUP catalogue 
is just a scaled-down version of the GALAXY catalogue, since 
GROUP is just a condensation of the GALAXY catalogue. In Fig. [3 
we plot the re-adjusted peculiar velocities (through Eq. ®) 
of two catalogues on the sky in the Galactic coordinate, so 
the Galactic plane region is empty. The red (light grey) points 
are moving away from us, and the blue (dark grey) points are 
moving towards us. The size of the points is proportional to the 
magnitude of the line-of-sight peculiar velocity. One can see 
that the velocity distribution is almost uniform across the full 
sky. 

’ The re-adjusted peculiar velocities of Cosmicflows-2 catalogue can 
be requested from the corresponding author. 



r [Mpc] 


Fig. 4. Covariance matrix of Vnir) for GALAXY catalogue. The .t- and 
y- axes mean the different radial bins, and the unit of colour bar is 
100 km s^f 


3. Methodology 

In this section, we first present three theoretical models of pair¬ 
wise velocity field, and then discuss how can we simulate the 
covariance matrix of the pairwise velocity estimator at different 
distance bins. Then we present the likelihood function. 
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3.1. Theoretical models 


Equation O is the definition of pairwise velocity field at each 
separate distance r. The approximate solution of the pairwise 
ve locity field through the pair conservation eq uation derived 
by Juszkiewicz et alJ (Il999h is given as (also see iFeldman et al.l 

l200l 


vnir) 

1(0 




f 


^{x)x^dx 




(5) 

(6) 


where a — 1.2 - 0.65y, and y = -(dln^/dlnr)|f=i. ^(r) can 
be measured directly from galaxy surveys two-point correla¬ 
tion statistics. Here we use the two-point correlation function 
of Point Source Catalogue re dshift (PSCz survey) calculated by 
iHamilton & Tegmartd (l2002h . This is because the PSCz survey 
measures the galaxy distance out to 120 fi^^Mpc and its sky cov¬ 
erage is close to t he Cosmicflows- 2 catalogue (see panels a, b, 
and d of fig. 1 in iMa et alJ (l2012h '). Thus the PSCz covers the 
similar volume of the local Uni verse as Cosmicflows-2. The ^(r) 
model derived by PSCz survey dHamilton & Tegmarkll2002h is 


^(0 = 




(7) 


where ri - 2.33 fi^'Mpc, r 2 - 3.51 fi^'Mpc, yi - 1.72, 72 = 
1.28, and erg is a free parameter. 

We name the above PSCz model as Model 1 in the following 
discussion. In addition, we also con sider two other simp lified 
models of the correlation function as (iFeldman et al.ll2()()^ 



where (7oTo)=(l-3, 4.76/!“'Mpc) and (1.8, T.fifi^’Mpc) for 
Models 2 and 3, respectively. 

We vary two parameters in the pairwise velocity field model. 
One is erg, which is the square root of the amplitude of the two- 
point correlation function (Eqs. © and ®). But because f(r) 
affects both the numerator and denominator of ^(r) (Eq. ®), erg 
will affect the shape of the pairwise velocity field function vi 2 (r). 
In addition, we also vary the total amplitude parameter as 
the pre-factor in Eq. ® . 

The three model predictions plotted in Fig.|5]are all negative 
in the regime of 0-30 Mpe, because of the gravitational potential 
effect. In addition, one can see that Models 1 and 2 predictions 
are close to each other, while Model 3 has slightly lower power 
on a large separation r ~ 30 Mpc. 


3.2. Covariance matrix 

We now calculate the covariance matri x for the pairwise ve - 
locity estimator. In the original work of iFerreira et alJ (Il999l) . 
the covariance matrix is calculated analytically by assuming a 
smoothed survey window function. However, this is only an ap¬ 
proximation since the real survey window can have complicated 
geometry. Here we calculate the covariance matrix of pairwise 
velocity field directly from a numerical simulation that automat¬ 
ically includes the real survey geometry. 

We simulated a large number of mock catalogues. For each 
simulation, we simulated A^data = 8315 and 5224 samples’ line- 
of-sight velocities for GALAXY and GROUP catalogues and then 


computed the pairwise velocities from Eq. ©. Then we calcu¬ 
lated the variance of the pairwise velocity field from these mock 
catalogues. We tested that Wim = 10^ is large enough for the 
results to converge. 

3.2.1. Covariance matrix of line-of-sight velocities 

As one can see from Eqs.(14)-(16) of iMa & Scott! (l2013h . the 
line-of-sight velocities are correlated. We therefore need to cal¬ 
culate the covariance matrix of these line-of-sight velocities be¬ 
fore simulation. For any two samples, their covariance matrix is 

Gij = H- Sijioi + aj) 

= {{rr\{r-;)){rj-\{rj))) + 6ij(aj+ 0-1), ( 9 ) 

where the first term is the cosmic variance term since the in¬ 
trinsic correlation between velocities at two different directions. 
The second bracket contains the nns measurement noise for each 
galaxy cr,- (i.e. quantity plotted in the upper panel of Fig.|2]l and 
the small scale and int rinsic dispersion (Jt which is estimated to 
be around 200 km s ' (iTurnbull et al.ll20l% . 

For the cosmic variance term, any two samples are corre¬ 
lated so the off-diagonal element of Gij (i + j) is non-zero. But 
for the measurement noise term, small scale, and intrinsic dis¬ 
persion, the two samples are not correlated, so they only con¬ 
tribute to the diagonal term in the Gij matrix. The first term is 
the real space velocity correlation function , which is related to 
the matter power s pectrum in Fourier space ( Watkins et al.l200^ 
lMa&Scottll2ni3h . 

((f, -v(r,))(f, -v(r,))) = ^4__- J dkPik)Fij{k), (10) 

where P(k) is the matter power spec trum that we output from 
public code cams (iLewis k al.ll2000h . The f{z) = din D/d In a 
is the growth rate function that characterizes how fast the 
structures grow at different epochs of the Universe. Since the 
Cosmicflows-2 samples peak at the redshift z ^ 0.0167 (lower 
panel of Fig. 13, we use the zero-redshift growth funct ion 
f{z = 0) = jWatkins et alJlTOO^ iMa & ScotlllMl in 

Eq. (fTOl) . In the future, if the survey probes deeper region of the 
space, one should use the corresponding growth function fiz) 
in Eq. flOll, so that the joint constraints on f(z)o-^ can be ob- 
tain ed, which constitutes a sensitive test of modify gravity mod - 

els (iHudson & Turnbuiill2012t[Planck Collaboration XIIli2015h . 

The Fij(k) is the window function 

— (f; ■ k) (f y ■ k) X exp(/kk ■ (r,- - r^)), (11) 

which can be calculated analytically (The Appendix in iMa et ^ 

IMTl) . 

Therefore, by calculating Gij matrix, we obtain a A^data x A^data 
covariance matrix for the line-of-sight velocities of mock galax¬ 
ies. We then followed the procedure in Appendix lAl to simulate 
a mock line-of-sight velocity catalogue. We did this repeatedly 
for Asim number of mock catalogues. 

3.2.2. Covariance matrix for pairwise velocities 

For each mock catalogue, we can plug them into Eq. (|2| to obtain 
viiir) for different bins, then one can obtain Asim numbers of 
mock V 12 velocity fields. Then, for each distance bin, one can 
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Fig. 5. vi 2 (r) for the GALAXY and GROUP catalogues with the data points 
calculated from Eq. (EJ and error bars from the square root of Eq. tnii. 
The black, red dashed, and blue dashed lines are for Models 1, 2, and 
3, respectively, by using Planck 2015 best-fitting cosmological param¬ 
eters. The measured Vi 2 (r) is separated into 14 bins, and the error bars 
between each bin are highly correlated. 


calculate the covariance matrix for A^bin as 

Cij = <vi2(r,)vi2(r;)) 

N ■ 

^Vsim 

Therefore, this C,j is a A^bin x M)in positive-dehnite symmetric 
matrix which is the covariance matrix for vi 2 (r). 

In Fig. m we plot the covariance matrix of the 
GALAXY catalogue. One can see that it is a positive and symmetric 
matrix, and the higher value of correlation exists on large separa¬ 
tion distances. This is because the cosmic variance term becomes 
more signihcant at a larger separation distance. In Fig.|5] we plot 
the square root of the diagonal value of covariance matrix as the 
error bars for GALAXY and the GROUP catalogue as an example 
for A^bin = 14. But here we remind the reader that the correlation 
between different data points is signihcant, so the value of the 
square-root cannot fully represent the total error budget. 

3.3. Likelihood 

We now have the models, measured peculiar velocity held data, 
and the covariance matrix ready. Our purpose is to ht the am¬ 
plitude of matter Huctuation parameter erg and the combined pa¬ 
rameter with the pairwise velocity data. We formulate a 
log-likelihood function as 

-\ogm = J]{v\2(n;0)-vUn))cr/ 

u 

X (v‘i2(o;«)-v?2(o))- (13) 

where the indexes “t” and “d” mean “theory” and “data” respec¬ 
tively, and 6 represents the parameters of interest. Maximizing 
this likelihood function will give us the estimate of the cosmo¬ 
logical parameter 6. 


4. Results 

In Fig. |5] the three models are plotted with Planck 2015 best- 
htting cosmological parameters, which seems to be consistent 
with both GALAXY and GROUP data sets. However, the error bars 
in Fig.|5]are highly correlated (Fig.lH, so we need to use the full 
covariance matrix to derive more quantitative results. 

We now present our hnal results for the likelihood anal¬ 
ysis. In Fig. |6] we plot the joint constraints on and erg 


from GALAXY and GROUP catalogues, for three theoretical mod¬ 
els. The black-diamond and orange-square marks are the best¬ 
fitting parameters from the WMAP nine-year results and Planck 
2015 results. One can see that the joint constraints from the 
GALAXY catalogue (upper panels) contain the best-htting point 
of the cosmic microwave background (CMB) survey at around 
Icr CL, which indicates that they are highly consistent with the 
ACDM model with WMAP or Planck cosmological parameters. 
For the GROUP catalogue (lower panels), although the best-fitting 
parameter values from CMB are not within the Icr CL of the joint 
constraints, they are consistent with pairwise velocity constraints 
at 2cr CL, which is n ot big enough to claim any discrepancy. In 
addition, as shown in IValentino et al.l (l20L5h . if one releasing 12 
cosmological parameters in Planck likelihood code, the crg-Q^ 
tension could be resolved. 

Comparing the six panels, we conclude that the current 
tightest constraint on cosmological parameters is from the 
GROUP catalogue with parameter values as 

= 0.102™, erg = 0.39!°;P (forModel 1). (14) 

Looking into the future, the potential for enhancing the con¬ 
straining power of the pairwise velocity held is embodied in the 
following aspects. 


1. Better modelling of the pairwise velocity from simulations 
of the l arge-scale str u cture. By developing a numerical sim¬ 
ulation, ISlosar et"^ (l2006l) investigated the distribution of 
velocities of pairs at a given separation taking bot h one-halo 
and tw o-halo contributions into account. Later on. lLam et al.l 
(1201 Ih studied how primordial non-Gaussianity affects the 
pairwise velocity probability density function by using an 
analytical model and the N-bod y simulations. More recently, 
iThompson & Nagaminel (1201 2l) have investigated the veloc¬ 
ity distribution of dark matter halo pairs in large N-body 
simulations (250/!“'Mpc-l A^^Gpe) and examined the pair¬ 
wise halo velocities with high velocity bullet cluster samples. 
These efforts at more accurate modelling of pairwise veloc¬ 
ity held will continue and lead to the closer description of 
the observed pairwise velocity data. 

2. Better calibration of the distance estimate and the reduc¬ 
tion of distance error. The most important aspect of enhanc¬ 
ing cosmological constraints on a pairwise velocity held 
is to i mprove the distance e stimate. Currently, the Tully- 
Fisher (iTullv & Fishei 19771) and FP (iDiorgovski & Davi^ 
ll987^l^mDbell et al.ll2014j) relations lead to the distance es¬ 
timate in the range of 25-30 per cent. For instance, for the 
FP distance estimator, the total scattering of distance r con¬ 
stitutes intrinsic scatters (~ 20 per cent), FP slope multiplied 
the observational error (~ 18 per cent) and the photomet¬ 
ric error (~ 3 per cent). The bulk part, i.e. intrinsic scatters, 
is due in part to the effect of stellar population age varia¬ 
tions on M/L, but the very large uncertainties on individual 
age estimates mean that it is very difficult to correct the ef¬ 
fect (ICampbell et al.lU014h . Therefore, much effort needs to 
be devoted to investigating those effects that strongly affect 
intrinsic scattering. 

3. Denser and broader sky coverage. These two factors in to¬ 
tal will lead to providing more samples for each individ¬ 
ual distance bin so that the nns noise level will be reduced 
as the sample size increases. In addition, since the Tully- 
Fisher survey needs to measure the HI line width, they are 
limited by the flux o f HI measurement (iMasters et al.ll200CT 
ICatinella et al.ll200^ . The Fundamental Plane survey mea¬ 
sures the galaxy’s surface brightness, so it is limited by the 
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Fig. 6. Joint constraints on Qf^^h-crg parameters with GALAXY and GROUP catalogues for three different models. Re d, blue, and green da shed 
lines are the 1, 2, 3 cr confidence levels, respectivel y. The best-fitting cosmological parameters from VTMAP 9-year (IHinshaw et alJl2013h and 
PlancklOlS results dPlanck Collaboration XTTlll2ni .ih are erg, /;) = {0.272, 0.82, 0.704}, and {0.309, 0.816, 0.677), respectively. The WMAP 
9-year values are taken as the VKMAP9+eCMB+BAOH-//o, and Planck 2015 parameter values are taken as Planck TT-fTE-t-EE-flowP-flensing-l-ext, 
where “lowP” means low-f" polarization data and “ext” means BAO-l-JLA-t-77o. 


sensitivity of the telescope (telescope size, exposure time, 
atmospheric interference, etc.). For these reasons, the faint 
galaxies are hard to measure, and the density in the survey 
volume is limited. As a rough estimate of the current state-of- 
the-art 6dF galaxy surveyj_the_sam£lingdens^ at zi /2 - 0.05 
is p = 4 X 10“^/j^Mpc“^ ( Campbell et ahllToi^ . 

4. Deeper redshift and distance surveys. The survey volume is 
an important factor in improving the dark energy and modi¬ 
fied gravity constraints. This is not only because the deeper 
survey provides more samples, but also because it provides 
estimates of the growth of the structure at different redshifts. 
As a result, for a much deeper survey, one can separate the 
pairwise velocity samples into different redshift bins and 
constrain the f{z)cr% value in Eq. (fTOl l. The redshift evolu¬ 
tion of /(z)crg in the regime of z = 0-1 constitutes a sensi¬ 
tive probe of the modified gravity and dark energy, w hich 
has been demonstrated in Hudson & TurnbulO d2012h and 
iPlanck Collaboration XIIll (|2015|) . 


5. Alterna tive velocity probe of t h e kine ti c Sun yaev-Zeldovich 
effect dSunvaev & Zeldovichl Il972i Il980ll . The kinetic 
Sunyaev-Zeldovich effect (kSZ) is the secondary anisotropy 
of the CMB produced by the inverse Compton scattering 
of the moving electrons. The temperature anisotropy of the 
CMB is therefore proportional to the line-of-sight peculiar 
velocities of the electrons, i.e. AT/T = -(crj/c)f dlne(r)(v- 
n). Thus, the accurate measurement of the kSZ effect 
from CMB observations and the precise modelling of in- 
tergalactic gas can, in principle, lead to determining the 
pairwise velocity field without suffering from systematics 
from optical sur vey. The pa i rwise momentum of the tem¬ 
perature field in iHand et ^ d2012l) gave the first detection 
of the kSZ effect by using the ACT data. More recently, 


IPlanck Collaboration XXXVl3 d2015h have re-measured the 
pairwise momentum field by using Planck 2015 data and 
cross-correlated the temperature field with the reconstructed 
peculiar velocity field. For the first time, they find the 
3-3.7cr CL detection of this correlation, which indicates 
that the gas is extended much beyond the virial radii 
of dark matter halo, consti tuting direct evidence of the 
low-density, diffuse baryons d Planck Coll aboration XXXVllI 
l2015l:^ernandez-Monteag^o et al.i2015h ! The effort of ac- 
curately measuring kSZ effect will carry on and will improve 
the peculiar velocity field measurement from the channel dif¬ 
ferent from optical surveys. 


5. Conclusion 

In this work, we have developed a Bayesian statistics method 
for using a galaxy pairwise velocity field to estimate cosmologi¬ 
cal paramete rs. We first reviewed t he pairwise velocity estimator 
developed in iFerreira et al.l (Il999h to calculate the relative mo¬ 
tion between galaxy pairs with a separation distance between 
0 and 30/r 'Mpe, and then we review three theoretical models 
of the mean pairwise ve locity field, especially the one used in 
lJuszkiewicz et akl (Il999l) . which was calibrated against N-body 
simulation. 

We then focused our effort on simulating mock galaxy cata¬ 
logues and calculated the log-likelihood function of the pairwise 
velocity field. We first calculated the covariance matrix between 
any two line-of-sight velocities in a given survey and simulated 
the 10^ mock catalogue according to this covariance matrix. 
Then we calculated the pairwise velocity field for each of these 
mock catalogues and calculated the variance of the vi 2 (r). In this 
way, we obtained a covariance matrix of the v\ 2 {r) field and used 
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this in the likelihood to estimate cosmological parameters 
and cTg. By using the Cosmicflows-2 catalogue, which consists 
of over 8 000 galaxies over the full sky out to 25 OOOkms^', we 
show that the results are consistent with the WMAP nine-year 
data and Planck 2015 data, and the error bars are currently too 
big to draw any concrete conclusion on modified gravity models 
or dynamical dark energjQ. We hope therefore that future data 
with more samples on the local Universe may help to pin down 
the sample variance and improve pairwise velocity field statis¬ 
tics. 
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Appendix A: Simulating a catalogue given its 
covariance matrix 

The Gjj is a symmetric, positive-definite matrix, where we can 
use the singular value decomposition (SVD) method to simulate 
the velocity vector v, (/ = 1, ...,Ndata)- The G matrix can be de¬ 
composed into 

G^U-S-U'^, (A.l) 

where U is the left matrix, and its transpose is the right ma¬ 
trix, and the S matrix is the diagonal matrix with Adata element 
that contains the eigenvalues of matrix G. Since G is a positive- 
definite matrix, the elements in S should be positive. Then we 
can simulate an array of Gaussian random variables v with zero 
mean and variance 5/ (i = 1 , 2 ,..., Adata) to satisfy 

(vv^) = S. (A.2) 

Finally, the real velocity vector we want to obtain is v = Uv, 
because 

(vv^> = {U ■ vv^ • U'^} 

= U ■ (vv^> ■ U'^ 

= U-S ■U'^ 

= G; (A.3) 

i.e. the covariance matrix of v is G. In this way, we simulate a 
line-of-sight velocity vector with dimension Adata, which satis¬ 
fies the desired covariance matrix. 
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